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We discuss some general aspects of acoustic black holes. We begin by describing the associated 
formalism with which acoustic black holes are established, then we show how to model arbitrary 
geometries by using a de Laval nozzle. It is argued that even though the Hawking temperature of 
these black holes is too low to be detected, acoustic black holes have interesting classical properties, 
some of which are outlined here, that should be explored. 
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I. INTRODUCTION 



Black holes are among the most fascinating objects in 
physics. The fact that they are pure objects, in the sense 
that they are made from spacetime itself, explains why 
they have taken such a special place in general relativity. 
There is a powerful and elegant mathematical machinery 
to describe them |l|, |2| , and their classical and quantum 
properties are well understood within the general relativ- 
ity framework. In this setting, the properties of isolated 
black holes have been thoroughly investigated. The much 
more complex processes that take part in the surround- 
ings of astrophysical black holes, the interaction of black 
holes with matter (accretion disks, magnetic fields, etc) 
Q or even with other black holes (for example, the prob- 
lem of the head-on collision of two black holes is now 
solved Q) are more or less well understood as well. On 
the semi-classical side, Hawking [BJ showed that when 
quantum effects are taken into account, black holes are 
not really black: they slowly evaporate by emitting an al- 
most thermal radiation. Hawking's prediction has been 
theoretically confirmed time and again in very different 
ways. The discovery of Hawking radiation uncovered a 
number of fundamental questions: among them the in- 
formation puzzle, the issue of the black hole final state, 
and so on. Some of these issues can be tackled only in a 
more fundamental theory, because classical general rela- 
tivity is not the ultimate theory of gravity since it does 
not embody the principles of quantum mechanics. A con- 
sistent theory of quantum gravity requires a modification 
of classical general relativity, and in the two alternative 
theories more in fashion nowadays, string theory and loop 
quantum gravity, black holes still occupy a special posi- 
tion. String theory's charm, for instance, derives in part 
from a couple of remarkable breakthroughs in connection 
with black hole physics, namely the entropy calculation 
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by a counting of micro states, and the computation of 
greybody factors |E S S @ • 

The progress in understanding black holes has been 
immense, over these last forty years since their concept 
was born, and they now play a central role in modern 
physics. Despite this, the lack of experimental tests has 
always been a drawback, for general relativists, and for 
people studying black holes in particular. An important 
step to make black holes more accessible (from an ex- 
perimental point of view) was given in 1981 by Unruh 
[lfj, who came up with the notion of analogue black 
holes. While not carrying information about Einstein's 
equations, the analogue black holes devised by Unruh do 
have a very important feature that defines black holes: 
the existence of an event horizon. The basic idea be- 
hind these analogue acoustic black holes is very simple: 
consider a fluid moving with a space-dependent velocity, 
for example water flowing throw a variable-section tube. 
Suppose the water flows in the direction where the tube 
gets narrower. Then the fluid velocity increases down- 
stream, and there will be a point where the fluid velocity 
exceeds the local sound velocity, in a certain frame. At 
this point, in that frame, we get the equivalent of an 
apparent horizon for sound waves. In fact, no (sonic) in- 
formation generated downstream of this point can ever 
reach upstream (for the velocity of any perturbation is 
always directed downstream, as a simple velocity addi- 
tion shows). This is the acoustic analogue of a black 
hole, or a dumb hole. These objects share more prop- 
erties with true, gravitational black holes, besides the 
existence of horizons: they display geodesies, wave ef- 
fects in their vicinity and, as we shall see they also emit 
Hawking radiation. Nevertheless they are not true black 
holes, because the acoustic metric satisfies the equations 
of fluid dynamics and not Einstein's equations. One usu- 
ally expresses this by saying that they are analogs of 
general relativity, because they provide an effective met- 
ric and so generate the basic kinematical background in 
which general relativity resides. They are not models 
for general relativity, because the metric is not dynam- 
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ical ly d ependent on something like Einstein's equations 
[Til [T^ . Following on Unruh's dumb hole proposal many 
different kinds of analogue black holes have been devised, 
based on condensed matter physics, slow light etcetera 
[Til IT2I Il3) . Analogue black holes have been the sub- 
ject of intense study because of the Hawking radiation 
they emit. In fact, it is now clear that the appearance 
of Hawking radiation does not depend on the dynamics 
of the Einstein equations, but only on their kinematical 
structure, and more specifically on the existence of an ap- 
parent horizon |l2l Il4j . The experimental verification of 
the Hawking effect is not easy though. Unfortunately as- 
trophysical black holes, having a Hawking temperature 
much smaller than the temperature of the cosmic mi- 
crowave background, accrete matter more efficiently than 
they evaporate. However, since Hawking radiation cru- 
cially depends on the existence of an apparent horizon, 
the analogues just described do emit Hawking radiation, 
and this was and still is the primary reason to study 
them. At present the Hawking temperatures associated 
to these analogues are too low to be detectable, but the 
situation is likely to change in the near future (see for 
instance 0). 

The importance of classical properties of analogue 
black holes have been somewhat underestimated. First, 
even though building (Hawking) very hot analogue black 
holes may be extremely difficult, building them any 
acoustic black hole is not. Thus we can easily have an 
acoustic black hole in almost any lab. What good are 
these black holes for? 

(i) They have an horizon. Hawking radiation is not the 
only interesting thing going on when an event horizon 
shows up! In particular, I would say that the "only in- 
going waves at the horizon" boundary condition would be 
interesting to see experimentally, with all its associated 
phenomena, some of which are mentioned below. 

(ii) Geodesies. This is a particularly interesting appli- 
cation. As will be shown one can easily mimic several 
geometries simply by varying the cross section of a do 
Laval nozzle. Thus we can observe the geodesies in dif- 
ferent spacetimes easily. 

(hi) Measuring absorption cross-sections. This would 
also be an interesting application of analogue black holes, 
to measure absorption cross-sections, glory effects, etc, 
and compare them with theoretical predictions |Tt| . 

(iv) Superradiance. This phenomenon, involving rotat- 
ing black holes was in the basis for the discovery of 
Hawking radiation [lg ■ To our knowledge this effect was 
never experimentally verified (not including Cherenkov 
radiation in this category |l9|). but it should not be very 
hard to rep roduce in the lab using acoustic black holes 

[Hil if. 

(v) Quasinormal modes. The resonance modes of black 
holes, called quasinormal modes (QNMs) are a very im- 
portant concept in any discussion involving gravitational 
radiation by black holes, the approach to equilibrium and 
black hole detection [23 0. The QNMs of black holes 
use the in-going waves at the horizon boundary condi- 



tion, and they usually have quite an interesting spectra 
of frequencies [2{|. The QNMs of some analogue black 
holes have already been computed [53, El H(j . 

(vi) Late-time tails. Black holes have no hair, and it 
is lost at late times as a power-law falloff [27| . Late-time 
tails can also be studied using black hole analogues [2(| • 

(vii) Analog black branes. It should be rather easy to 
implement other analogue black objects, such as black 
branes and strings, which could deepen our understand- 
ing about these geometries. 

(viii) Interaction of black holes with electric and mag- 
netic fields. On a more speculative vein, it is even pos- 
sible in principle to simulate in the lab the interaction 
of astrophysical black holes with matter and with elec- 
tromagnetic fields. It is even possible that one might be 
able to study effects such as the Blandford-Zjanek pro- 
cess 0, • This would be a tremendous motivation to 
use and explore analogue black holes. Some steps along 
this direction, although not directly connected to ana- 
logue black holes, were given in [29) . 

These are just classical aspects of black holes, but even 
these must be mastered before embarking on experimen- 
tal Hawking radiation detection. Not only must one con- 
trol what happens in the experimental situation, but the 
understanding of classical phenomena may bring clues on 
how to favor the probabilities to detect Hawking radia- 
tion. It is also worth stressing that some purely classical 
phenomena shed light on quantum aspects of (analogue 
and general-relativistic) black hole physics. For exam- 
ple, positive and negative norm mixing at the horizon 
leads to non-trivial Bogoliubov coefficients in the calcu- 
lations of Hawking radiation |30| ; superradiant instabili- 
ties of the Kerr metric are related to the quantum process 
of Schwinger pair production 0, [32, H3, 03 j an d more 
speculatively (classical) highly damped black hole oscil- 
lations could be related to area quantization [35) (this 
possibility was discarded in [2(j because of the failure 
to satisfy the laws of black hole mechanics j3|| , but the 
situation may change [37) 1- 

Most of what has been said applies equally well to other 
types of analogue black holes (see for instance [T^. I3S)l 
but for simplicity we shall here deal only with acoustic 
black holes. Some aspects of acoustic black holes will be 
explored: we'll explain how to generate a large class of 
acoustic black holes by using a de Laval nozzle with a 
variable cross section profile. We will see how to make a 
simple black brane and study its stability properties. We 
will then explain wh y in certain situation the analogue 
branes are unstable [39) . We will make a small review 
of what has been done so far concerning classical aspects 
of wave propagation in acoustic black holes, taking the 
(2 + l)-dimensional rotating black hole as a case study. 



II. EFFECTIVE ACOUSTIC GEOMETRY 

This section will be as self contained as possible, be- 
cause we want to make explicit the assumptions that go 
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with the usual derivation of the acoustic metric. How- 
ever, this derivation can be found in the monograph by 
Matt Visser . Let us start with the equations of fluid 
dynamics, and try to arrange them in such a way that an 
effective metric stands out naturally. The fundamental 
equations of fluid dynamics |40l Ell E^ are the equation 
of continuity 



d t + v.(H = o, 

and Euler's equation 



(1) 



p-r = p[dtv + (v.V)v] — —Vp + F , (2) 
at 

where F are for the moment all the external forces act- 
ing on the fluid. Hereafter we shall make the follow- 
ing assumptions: (i) the external forces are all gradient- 
derived, or F — — pV<I>. Thus we are neglecting vis- 
cosity terms in Navier-Stokes equation, (ii) the fluid is 
locally irrotational, and introduce the velocity potential 
ip, v = — V^; and (iii) the fluid is barotropic, i.e., the 
density p is a function of pressure p only. In this case, 
we can define 



h(p) 



V/i = 



dp' 
P(P') 



Euler's equation can be written as 



1 



d t i> + h + ^{Vi>f + $ = 0. 



(3) 



(4) 



(5) 



To study sound waves, we follow the usual procedure 
and linearize the continuity and Euler's equations around 
some background flow, by setting p = po + ^Pi > P = Po + 
epi , V — V'o + £"01 > an d discarding all terms of order e 2 
or higher. The external potential $ is taken as constant. 
The continuity equation yields 

ftpo + V.(po«o) = 0, (6) 
d t pi + V.(pi« + PqVi) = 0. (7) 

Linearizing the enthalpy we get h(pg + epi) ~ h(po) + 
e^r = ho + e— . Inserting this in Euler's equation 

a Pp—PQ PO 

one gets 

1 



9^o + ^o + ^(V^o) 2 + * = 0, 



(8) 

Pi = Po{d t ipi + v .VV>i) = • (9) 
On the other hand, since the fluid is barotropic we have 

dp 



Pi = ttPi , 
op 

and using © this is the same as 
dp 

Pi = -£-pa{dti>i + vo-V^i) ■ 
op 



(10) 



(11) 



Finally, substituting this into J7J we get 

o- - ( dp 



d t \-g^Po{d t i/>i + Vo-V^i) 
dp 

V. [ poV-01 - -7^PoVo(dtipi + v .VV>i) 



(12) 
(13) 



It can now easily be shown \H\ that this equation can also 
be obtained from the usual curved space Klein-Gordon 
equation 



with the effective metric g^ v given by 



1 

Poc 



(14) 



(15) 



Neither of the background quantities is assumed constant 
through the flow, and so they are in general dependent 
on the coordinates along the flow. Here we have used 
the definition of the local sound speed c~ 2 = We 
can see that the propagation of sound waves in a fluid is 
equivalent to the propagation of a scalar field in a generic 
curved spacetime described by (|15|l . or in covariant form 
by 



= p0 

c 



(16) 



This means that all properties of wave propagation in 
curved space hold also for the propagation of sound 
waves. The power of this effective geometry should be 
clear: first, by changing the background flow, we change 
the effective acoustic metric. Second, since this geome- 
try clearly has an apparent horizon at the point where 
Vq = c, and since the existence of an apparent horizon 
implies Hawking radiation, then there should be Hawk- 
ing radiation in this geometry, which takes the form of 
phonons (To) . The Hawking temperature can be com- 
puted to yield [TT| : 

h djc-vjj 

kT H = ^ , 17 

2tt on 

where v± is the component of the fluid velocity normal 
to the horizon, and n is the unit vector normal to the 
horizon. This can also be written as 



Tjj = 1.2 x 10 _9 Km 



c 




1 d(c — v±) 


llOOOms- 1 - 




c dn 



(18) 



This is, for all practical purposes a number too low to 
be detected, and it is even more so if one notices that 
one has to deal with the ambient noise. Despite the fact 
that one cannot observe Hawking radiation, one can still 
measure classical aspects of black holes. So we now turn 
to this, but first we explain how we can mimic several 
geometries by using a de Laval nozzle. 
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v < c. 




FIG. 1: A sketch of a de Laval nozzle, used to make a smooth 
transition from subsonic to supersonic flow. The velocity of 
the fluid v(x) equals the local velocity of sound c(x) at the 
throat of the nozzle, x = xt- The cross-section at this point 
is denoted by At- 



III. SHAPING THE NOZZLE 

A. The de Laval nozzle 

A de Laval nozzle is a device which can be used to 
accelerate a fluid up to supersonic velocities. They were 
first used in steam turbines, but they find many applica- 
tions in rocket engines, nozzles in supersonic wind tun- 
nels, etc. It consists of a converging pipe, where the fluid 
is accelerated, followed by a throat which is the narrow- 
est part of the tube and where the flow undergoes a sonic 
transition, and finally a diverging pipe where the fluid 
continues to accelerate. It is sketched in Fig. ^ 

Consider now a steady, isentropic flow through the noz- 
zle, which has a varying cross-section A(x), where x is 
the arc length along a streamline. Logarithmic differen- 
tiation of the continuity equation 



yields 



dm 

pvA — —— = const . 
at 



1 dv 1 dA I dp 

- v T x + Ad X - + ~d=°- 



(19) 



(20) 



For isentropic flow p = p(p) and the definition for the 

immediately gives 



speed of sound 

„2 



dp 

dp 



constcntropy 



& . Using this in we obtain 



1 dv 

v dx 



1 dA 

A dx 



1 dp 

c 2 p dx 



0. 



(21) 



Using the component of Eulcr's equation along the 
streamline 



pv 



dv 
dj- 



dp 

dx ' 



(22) 



and combining it with equation (|21|l we have finally 



:(! 



v . dv 
c 2 dx 



1 dA 
A dx 



(23) 



According to this differential equation, when the flow is 
subsonic (v < c), dv / dx and dA/dx have opposite signs. 



So narrowing the pipe will make the gas flow faster, 
which is what we expect from common experience. In 
fact, for very small v, equation (|23[) can be written as 
dv jv — —dA/A and thus vA is a constant, a well known 
result for incompressible fluids. The situation is opposite 
for v > c, when dv/dx and dA/dx have the same sign. 
This means that a region of increasing cross- section will 
accelerate the flow. The nozzle equation (|23H shows that 
the transonic flow through the nozzle must reach v = c 
at the throat where dA/dx = 0. This is a necessary con- 
dition, but not a sufficient one. Whether the actual flow 
will be transonic depends of course on the lower bound- 
ary condition, i.e., the velocity vq at the entrance of the 
nozzle. For a given profile A(x) of the nozzle v$ has to 
have exactly the right value: if vq is too small the flow 
will remain subsonic everywhere. If vq is too large the 
velocity will reach v = c upstream from the throat, at 
x < xt, where dA/dx < and so dv/dx — > — oo at that 
point. The flow will stagnate which results in a shock 
between the flow and the low velocity region. So in this 
case there is no smooth transonic flow. 

I shall from now on assume that the boundary con- 
ditions are such that there is a smooth transition from 
sub to supersonic flow, at the throat located at x = xt- 
Given a nozzle profile A(x) and an equation of state, then 
Eq. 12MII together with H16|) describes completely, for 
x < xt, an acoustic black hole (I note that the equation 
of state p = p{p) will allow, from Eq. (|20() . to have c as 
a function of v). Conversely, given an equation of state, 
a whole family of acoustic black holes can be obtained, 
by simply varying A(x). For example, for a perfect gas, 



ilyvi 

we have |4(| c 2 = c 2 , — (7 — l)/2w 2 , where cq is the speed 
of sound at the location corresponding to v = 0, and 7 
is the ration of specific heats. Now that we have c as 
function of v, Eq. (|23[) allows one to specify the velocity 
profile and therefore the full metric (|16|) as a function of 
A(x). 



B. Acoustic black holes by different nozzle 
configurations 

Let us now suppose that the dependence of c on x is 
small, i.e., let us assume that dc/dx is not very large 
(as happens for example for perfect gases) , and therefore 
that c ~ const. With the assumption of constant sound 
velocity, equation (|2*5|) can be solved yielding 



= « 2 /(2c 2 ) 



A{x)C . 



(24) 



where the constant C is found by applying the condition 
at the throat v(xt) = c. This gives 



1 



y/(2c 2 )-l/2 



A(x)/A(x T ) 



Let us choose the following generic form for A: 
1 



A(x) 



D [A(x T )f(x)Y/2-l/2 



(25) 



(26) 
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For this to be a consistent solution we must have 
dA/dx — at x = Xt, which results in the constraint 



f'(x T )e lA( - XT)nxT)]2/2 - 1/2 



A(x T ) 2 --—2 =0. 

(27) 

This constraint, together with the definition (|26|) can be 
satisfied if one chooses 



Equation i|25|) is then trivially solved by 
- = I(x)A{x T ) 



(28) 



(29) 



So we conclude that in order to mimic some metric, we 
have to be able to simulate the correct background flow, 
as indicated by Eq. I|16l) . Now, to mimic the background 
velocity, one only has to build a de Laval nozzle according 
to H29JI . and we have our problem solved. One can also 
look at the Hawking radiation of acoustic black holes in 
nozzles. This was done recently by Barcelo, Liberati and 
Visser f43|. 



IV. CLASSICAL WAVE PHENOMENA NEAR 
ACOUSTIC BLACK HOLES 

Some classical aspects of wave propa gati on in acoustic 
black holes have been explored in [2(J, |21| . Here I will 
summarize some of their results and also comment on 
absorption cross-sections, focusing always on the (2 + 1)- 
dimensional rotating acoustic black hole [Tlj , which I now 
describe. 

A simple rotating acoustic black hole geometry was 
presented in [Tl| modeling a "draining bathtub", ideal- 
ized as a (2+l)-dimensional flow with a sink at the origin. 
I will show that a simple generalization of this geometry 
can mimic a black brane, and I will show why, despite 
recent claims |39j, it is not unstable. 

Consider a fluid having (background) density p. As- 
sume the fluid to be locally irrotational (vorticity free), 
barotropic and inviscid. From the equation of continu- 
ity, the radial component of the fluid velocity satisfies 
pv r ~ 1/r. frrotationality implies that the tangential 
component of the velocity satisfies v e ~ 1/r. By con- 
servation of angular momentum we have pv e ~ 1/r, so 
that the background density of the fluid p is constant. In 
turn, this means that the background pressure p and the 
speed of sound c are constants. The acoustic metric de- 
scribing the propagation of sound waves in this "draining 
bathtub" fluid flow is ED: 



ds 2 



A 2 + B' 2 



2Bdcj)dt + dr 



- r 



dt 2 

2 ji2 



2A 



-drdt 



Here A and B are arbitrary real positive constants related 
to the radial and angular components of the background 
fluid velocity: 



-At + B6 



(31) 



In the non-rotating limit B = the metric ()30|1 reduces 
to a standard Painleve-Gullstrand-Lemaitre type metric 
.45]. The acoustic event horizon is located at rjj = A/c, 
and the ergosphere forms at tes = {A 2 + B 2 ) 1 / 2 jc. 

Some physical properties of our "draining bathtub" 
metric are more apparent if we cast the metric in a Kerr- 
like form performing the following coordinate transfor- 
mation (where again we correct some typos in 



dt^dt- 



Ar 
r 2 c 2 — A 2 



BA 



r(r 2 c 2 — A 2 ) 



dr . 



(32) 



Then the effective metric takes the form 



ds 2 = -11- 



1 - 



A 2 



dr 2 - 2Bd4>dt + r 2 d4> 2 . (33) 



(30) 



Notice an important difference between this acoustic 
metric and the Kerr metric: in the (t, t) component of 
the metric Q33JI the parameters A and B appear as a sum 
of squares. This means that, at least in principle, there 
is no upper bound for the rotational parameter B in the 
acoustic black hole metric, contrary to what happens in 
the Kerr geometry. 



A. QNMs 

Black holes, like so many other objects, have character- 
istic oscillation or ringin g m odes, which are called quasi- 
normal modes (QNMs) [23, |2J, |25|, the associated fre- 
quencies being termed QN frequencies, or wqn- The QN 
frequencies of the (2 + l)-rotating acoustic black hole, 
described by the metric l|3*3*jl were recently computed in 
pol |20| , and so were the QNMs of the canonical acous- 
tic black hole. The numerical results, consistent with a 
WKB analysis, are shown in Figs. 12151 

(i) m > 0: In Fig. 12131 we show results pertaining to 
perturbations having positive m, i.e., co-rotating waves. 
In Fig. |2 we show the real part of the QN frequencies 
for m = 1 modes as a function of the black hole rotation. 
Higher m modes follow a similar pattern. One can see 
from this plot that for low black hole rotation parameter 
B the different overtones are clearly distinguished, but 
that as the rotation increases they tend to cluster and 
behave very similarly. For very large rotation B, all the 
overtones behave in the same manner, and in this high 
rotation regime the real part of the QN frequency scales 
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FIG. 2: The real part of the QN frequency as a function of 
the rotation parameter B/A, for several overtones of a m = 1 
mode. Here, th = A/c is the horizon radius. Note how all 
the several lowest overtones "coalesce" in the high rotation 
regime, all growing linearly with B/A. 




FIG. 3: The imaginary part of the QN frequency as a function 
of the rotation parameter B/A, for several overtones of a m — 
1 mode. It is clear from this plot that the imaginary part of 
the QN frequencies of m > modes is very insensitive to the 
rotation of the black hole. 



linearly with the rotation. Indeed we find that the slope 
is also proportional to m so that 



Re[aj QN ] 



mBc 2 
~A2~ 



D 



oo . 



for to > (34) 



We notice that this behavior was already present in the 
WKB investigation in [2(|. In Fig. we show the imag- 
inary part of the QN frequencies as a function of the 
rotation parameter, for m = 1. Different overtones have 
different imaginary parts. Note also that for high B the 
real part of the modes coalesce whereas the imaginary 
part does not. The magnitude of Im[wQjv] increases with 
B, which was observed also in the WKB approach [2fj| . 
Thus, as the rotation increases the perturbation dies off 



FIG. 4: The real part of the QN frequency as a function of the 
rotation parameter B/A, for several overtones of a m = — 1 
mode. Notice that for each overtone number n there is a 
critical rotation at which the mode crosses the axis, i.e., there 
is a critical rotation B/A at which the real part of the QN 
frequency is zero. Higher overtones cross the axis at a slower 
rotation. We have not been able to follow the mode beyond 
this point. 
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FIG. 5: The imaginary part of the QN frequency as a function 
of the rotation parameter B/A, for several overtones of a m = 
— 1 mode. We have not been able to follow the modes beyond 
a certain critical point (defined as the rotation B/A for which 
the real part of the QN frequency is zero) . Nevertheless, an 
judging by the modes we did manage to follow, namely the 
fundamental mode, it seems that Im[ojQjv] never crosses the 
axis, i.e., it is always negative, and therefore the mode is 
stable. 



quicker. This also means that the black hole is stable 
against m > perturbations, because the imaginary part 
is always negative. 

(ii) to < 0: In Figs. 1 11)1 we show results concerning 
perturbations having negative m, i.e., counter- rotating 
waves. The behavior of the QN frequencies for m < 
is drastically different from the to > perturbations. In 
Fig. 0]we plot the dependence of Re[wQjv] as a function 
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of the rotation of the black hole B. As B increases the 
magnitude of the real part of the QN frequency decreases. 

The oscillation frequencies for the fundamental modes, 
labeled by n = 0, indeed get close to the horizontal axis as 
B goes to infinity. However, we haven't been able to track 
some overtone modes with negative m for very high rota- 
tion since, as can be seen in Fig. 0] the real part of these 
modes eventually change sign. It is extremely difficult, 
using the method employed here, to compute modes hav- 
ing Re[wQjv] ~ 0. Nevertheless, supposing that (as the 
numerical studies for the fundamental modes indicate) 
the QN frequencies asymptote to zero for very large B, a 
WKB analysis shows that loqn ~ k/B, where k is some 
m-dependent constant. The imaginary part of the QN 
frequencies behaves in a similar manner, as seen in Fig. 

El 

(iii) m = 0: For circularly symmetric (m = 0) modes, 
our numerical method shows no sign of convergence. For 
in = 0, the wave equation can be written in the simpler 
form 



*,f.f. +(^ 2 -^)* = 0, 



where 



V 



f» _ i 



1 5 



(35) 



(36) 



The potential V is not positive definite, and this pre- 
cludes also a simple stability proof. 

To have a better physical understanding of this data, 
consider a m = 1 mode, for which the lowest mode (this 
is the mode that controls the ringing phase) is approx- 
imately loqn ~ (0.4 — 0.33i)^-, with m the horizon 
radius. If one builds an acoustic black hole by mak- 
ing a Th = 1 mm hole in a tub with water, then this 
black hole should have a characteristic ringing frequency 
of lo ~ 4x 10 5 s _1 , and a typical damping timescale given 
by r = t^tt ~ 3 x lCT 6 s. 



B. Late-time tails 

The existence of late-time tails in black hole spacetimes 
is well established, both analytically and numerically, in 
linearized perturbations and even in a non-linear evolu- 
tion, for massless or massive fields [27j- This is a prob- 
lem of more than academic interest: one knows that a 
black hole radiates away everything that it can, by the 
so called no hair theorem (see |46( for a nice review), but 
how does this hair loss proceed dynamically? A more or 
less complete picture is now available. The study of a 
fairly general class of initial data evolution shows that 
the signal can roughly be divided in three parts: (i) the 
first part is the prompt response, at very early times, and 
the form depends strongly on the initial conditions. This 
is the most intuitive phase, being the obvious counter- 
part of the light cone propagation, (ii) at intermediate 



times the signal is dominated by an exponentially decay- 
ing ringing phase, and its form depends entirely on the 
black hole characteristics, through its associated quasi- 
normal modes [23, [2^, [47| • (iii) a late-time tail, usually 
a power law falloff of the field. This power law seems to 
be highly independent of the initial data, and seems to 
persist even if there is no black hole horizon. In fact it 
depends only on the asymptotic far region. 

It is not generally appreciated that there is another 
case in which wave propagation develops tails: wave 
propagation in odd dimensional flat spacetimes. In fact, 
the Green's function in a D-dimensional spacetime (see 
Cardoso et al in [2]j and also 0> 03 ) have a completely 
different structure depending on whether D is even or 
odd. For even D it still has support only on the light 
cone, but for odd D the support of the Green's function 
extends to the interior of the light cone, and leads to the 
appearance of tails. 

Analogue black holes also exhibit late-time tails, shed- 
ding their hair in a power-law falloff manner. As an ele- 
gant application of the wave tail formalism developed by 
Ching et al , it was found in [2(j that any perturbation 
in the vicinity of the (2 + l)-dimensional analogue black 
hole described by ffify eventually decays as a power-law 
falloff of the form 



H ~ £-( 2 "»+l) 



(37) 



On the other hand, this s precisely the tails that appear 
in any (2 + l)-dimensional flat spacetime (see Cardoso et 
al in (23). We thus have a consistent and elegant result. 



C. Superradiant amplification of phonons 

Rotating black holes can superradiate, in the sense 
that in a scattering experiment the scattered wave has 
a larger amplitude (the frequency is the same, this is 
not a Doppler effect has explained in [l^ an d references 
therein) than the incident wave. Superradiance is a gen- 
eral phenomenon in physics. Inertial motion superradi- 
ance has long been known |5Cf . and refers to the pos- 
sibility that a (possibly electrically neutral) object en- 
dowed with internal structure, moving uniformly through 
a medium, may emit photons even when it starts off in its 
ground state. Some examples of inertial motion superra- 
diance include the Cherenkov effect, the Landau criterion 
for disappearance of superfluidity, and Mach shocks for 
solid objects traveling through a fluid (cf. for a dis- 
cussion). Non- inertial rotational motion also produces 
superradiance. This was discovered by Zel'dovich [l7| . 
who pointed out that a cylinder made of absorbing ma- 
terial and rotating around its axis with frequency £1 can 
amplify modes of scalar or electromagnetic radiation of 
frequency to, provided the condition 



lo < mfl 



(38) 



(where m is the azimuthal quantum number with respect 
to the axis of rotation) is satisfied. Zel'dovich realized 
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that, accounting for quantum effects, the rotating object 
should emit spontaneously in this superradiant regime. 
He then suggested that a Kerr black hole whose angular 
velocity at the horizon is f2 will show both amplification 
and spontaneous emission when the condition l|38|) for 
superradiance is satisfied. This suggestion was put on 
firmer ground by a substantial body of work |5lj . In par- 
ticular, it became clear that (even at the purely classical 
level) superradiance is required to satisfy Hawking's area 
theorem |52ll53| . 

Superradiance is essentially related to the presence of 
an ergosphere, allowing the extraction of rotational en- 
ergy from a black hole through a wave equivalent of the 
Penrose process [54j- Under certain conditions, superra- 
diance can be used to induce instabilities in Kerr black 
holes HjJ. Indeed, all spacetimes admitting an ergo- 
sphere and no horizon are unstable due to rotational 
superradiance. This was shown rigorously in but 
the growth rate of the instability is too slow to observe 
it in an astrophysical context pa, E3 ■ Kerr black holes 
are stable, but if enclosed by a reflectin g m irror they can 
become unstable due to superradiance |52t |58| ; 

The possibility to observe rotational superradiance in 
analogue black holes was considered by Schutzhold and 
Unruh [3 8 |. and more extensively by Basak and Majum- 
dar [53, |6(j, who computed analytically the reflection 
coefficients in the low frequency limit cuA/c 2 <§; 1. In 
particular, the authors of [3^1 showed that the ergore- 
gion instability in gravity wave analogues is related to 
the existence of an "energy function" [their Eq. (68)] 
that is not positive definite inside the ergosphere. In the 
context of analogues, inertial superradiance based on su- 
perfluid 3 He has been studied by Jacobson and Volovik 
|61j . The explicit numerical calculation of reflection co- 
efficients for the (2 + l)-rotating acoustic black hole was 
done in |20j . in the superradiant regime. 

Results of the numerical integrations for the draining 
bathtub metric are shown in Fig. |S| Panels on the left 
show the reflection coefficient |i? w i| 2 for m = 1, and pan- 
els on the right show |i?^2| 2 for m — 2, for selected values 
of the black hole rotation B. Panels on top show that, 
as expected, in the superradiant regime < uj < mB the 
reflection coefficient |i? wm | 2 > 1. Furthermore, as one in- 
creases B the reflection coefficient increases, and for fixed 
B, the reflection coefficient |i? wrra | 2 attains a maximum at 
lo ~ mB, after which it decays exponentially as a func- 
tion of ui outside the superradiant interval. This is very 
similar to what happens when one deals with massless 
fields in the vicinities of rotating Kerr black holes [52| . 
In particular, from the close-up view in the middle pan- 
els we see that, for B = 1, the maximum amplification is 
21.2 % (to = 1) and 4.7 % (to = 2). 

As a final remark, and as we have anticipated, an im- 
portant difference between the acoustic black hole metric 
and the Kerr metric is that in the present case there is 
no mathematical upper limit on the black hole's rota- 
tional velocity B. In the bottom panels we show that, 
considering values of B > 1, we can indeed have larger 
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FIG. 6: Reflection coefficient |i?u m | as a function of iv for 
m — 1. Each curve corresponds to a different value of B, as 
indicated. The top panels show that the reflection coefficient 
decays exponentially at the critical frequency for superradi- 
ance, ujsr = mB. The middle panels show a close-up view in 
the superradiant regime for B < 1: at B = f the maximum 
amplification is 21.2 %. The bottom panels show that su- 
perradiant amplification can become much more efficient for 
values of the rotation parameter B > 1. 



amplification factors for acoustic black holes. 

Summarizing: if we are clever enough to build in the 
lab an acoustic black hole that spins very rapidly, rota- 
tional superradiance can be particularly efficient in ana- 
logues. This is an important result, considering that the 
detection of rotational superradiance in the lab is by no 
means an easy task, as originally predicted by Zel'dovich 
|17| and confirmed by recent reconsideration of the prob- 
lem |Tg| . Of course, in any real- world experiment the 
maximum rotational parameter will be limited. At the 
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mathematical level, the equations describing sound prop- 
agation (which are written assuming the hydrodynamic 
approximation) will eventually break down. Physically, 
if the angular component of the velocity v s becomes very 
large the dispersion relation for the fluid will change, in- 
validating the assumptions under which we have derived 
our acoustic metric [2(J Hg • 

The superradiant phenomena we have described are 
purely classical in nature. However, an interesting sug- 
gestion to observe quantum effects in acoustic superradi- 
ance was put forward in |60| . To write down our acoustic 
metric we required the flow to be irrotational and non- 
viscous. As a natural choice, we could use a fluid which is 
well known to possess precisely these properties: super- 
fluid Hell. In this case the presence of vortices with quan- 
tized angular momenta may lead to a quantized energy 
flux. The heuristic argument presented in [6(j goes as fol- 
lows. Let us imagine that our black hole is a vortex with 
a sink at the center. In the quantum theory of Hell the 

wavefunction is of the form '5 = exp i ^ . (j)(fj ; )$ gTO und j 

where fj is the position of the j-th particle of Hell. The 
velocity at any point is given by the gradient of the phase 
at that point, v = V(/>, so that (roughly speaking) the 
velocity potential can be identified with the phase of the 
wavefunction. This phase will be singular at the sink 
r = 0. Continuity of the phase around a circle surround- 
ing the sink requires that the change of the wavefunction 
satisfies Acf) — 2irB. For the wavefunction to be single 
valued, B (that is, the black hole's angular velocity at the 
horizon) must be the integer multiple of some minimum 
value AB, i.e., B = nAB. Then the angular momentum 
of the acoustic black hole would be forced to change in in- 
teger multiples of AB. Correspondingly, the spectrum of 
the reflection coefficients may be given by equally-spaced 
peaks with different strengths. This discrete amplifica- 
tion could enhance chances of observing superradiance 
in acoustic black holes, and rule out (or provide empir- 
ical support to) some of the many competing heuristic 
approaches to black hole quantization. 



D. Acoustic black branes and superradiant 
instabilities 



Now, we are free to add an extra dimension and inter- 
preting the result as the superposition of a vortex fila- 
ment and a line sink ^lj . We get therefore the following 
line element 



ds 1 = - 1 - 



A 2 + B 2 
c 2 r 2 
2 a 12 



c 2 dP 



A 2 



dr 2 



2Bd<pdt + r z d(j) z + dz 



(39) 



This describes an analogue black brane, and compactifi- 
cation of the transverse direction z can be accomplished, 
in practice by using a "tamper" around the flow. 

The propagation of a sound wave in a barotropic invis- 
cid fluid with irrotational flow is described by the Klein- 



Gordon equation V M V A ' V I' = for a massless field ^ in 
a Lorentzian acoustic geometry, which in our case takes 
the form (|39|l . In our acoustic geometry we can separate 
variables by the substitution 



1 



(40) 



and we get the wave equation 



V = f ro 



Bm N 


2 


r 2 , 


) -\ 






~ 1) 


+ 4^ 



(41) 



(42) 



To arrive at (|42|) we have already performed the following 
re-scaling: f = rA/c, u) = uiA/c 2 , B = B/A, fi = [icjA. 
The re-scaling effectively sets A — c — 1 in the original 
wave equation, and picks units such that the acoustic 
horizon fjj = 1. The quantity / = (1 — -^4^), and the 
tortoise coordinate is defined by the condition 



dr* 
dr 



1 
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Explicitly, 



A. 

2c l ° g 



-A 



cr - 



A 



(43) 



(44) 



It is known that the Kerr geometry, or any rotating 
(absorbing) body displays superradiance 0,112 This 
means that in a scattering experiment of a wave with 
frequency u> < mfl the scattered wave will have a larger 
amplitude than the incident wave, the excess energy be- 
ing withdrawn from the object's rotational energy. Here 
n is the horizon's angular velocity and m is an azimuthal 
quantum number. Now suppose that one encloses the 
rotating black hole in a spherical mirror. Any initial per- 
turbation will get successively amplified near the black 
hole event horizon and reflected back at the mirror, thus 
creating an instability. This is the black hole bomb, as 
devised in |52j | and recently improved in [581 ]. This in- 
stability is caused by the mirror, which is an artificial 
wall, but one can devise "natural mirrors" if one consid- 
ers massive fields. Imagine a wavepacket of the massive 
field in a distant circular orbit. The gravitational force 
binds the field and keeps it from escaping or radiating 
away to infinity. But at the event horizon some of the 
field goes down the black hole, and if the frequency of the 
field is in the superradiant region then the field is ampli- 
fied. Hence the field is amplified at the event horizon 
while being bound away from infinity. Yet another way 
to understand this, is to think in terms of wave propaga- 
tion in an effective potential. If the effective potential has 
a well, then waves get "trapped" in the well and amplified 
by superradiance, thus triggering an instability. In the 
case of massive fields on a (four-dimensional) Kerr back- 
ground, the effective potential indeed has a well, as we 
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FIG. 7: A typical form for the effective potential in the Kerr 
geometry, here shown for I = m = 1 modes. We have set 
the mass of the black hole M = 1, so the rotation parameter 
a varies between (Schwarzschild limit) and 1/2 (extremal 
limit). Here we plot the effective potential for the near ex- 
treme situation, a ~ 0.5 and for fx — 0.7 and w = 0.6878. 



show in Figure Consequently, the massive field grows 
exponentially and is unstable (see 0, 13^, H3, HH for 
explicit examples). With this in mind, we would expect 
that black strings and branes of the form l|39|) for which 
there is a bound state will be unstable; here the trans- 
verse direction z works as an effective mass for the sound 
wave. To get a bound state, one necessary condition is 
that the derivative of the potential is positive, at asymp- 
totic large radial distances (see 39] for more details). 
Now, near infinity, we get 



Veff 

which leads to 



V'tt 



Am 2 — 1 — p 2 + 8uimB 
4^2 



4/i 2 + 1 - 4m 2 - 8u>mB 
2r~ 3 



(45) 



(46) 



Now, this can be positive, thus an instability can be trig- 
gered. 

For the sake of generality, let us drop the constant 
p requirement, which means not assuming conservation 
of angular momentum (this can be achieved by having 
external torques), then we can show that the effective 
metric is 



ds 2 = - 1 - 



A 2 /p 2 + B 2 
c 2 r 



1 - 



A 



> x -1 



pi c 2 r 2 



dr 2 - 2Bd4>dt + r 2 d^ 2 + dz 2 . 

(47) 



Here, A , B are again constants but they have different 
dimensions. Separating variables by the substitution 



(48) 



implies that 'F(r) obeys the wave equation (just insert 
the ansatz in Klein-Gordon's equation) 



eP* / Bm ^ 

+ ( {LJ -— ) 



V * = 0. 



(49) 



Here 



V = / U 2 c + 



Am 2 c — c 



c 

2^ 



/ 



A 2 



Ar 2 



5 + 2r( C - + 2 J±* 
c p ' 



^Acr 4 p 2 

and the tortoise coordinate is defined as 



dr 
dr* 



c(l 



A 2 



) = f- 



(50) 



(51) 



Notice that for constant p , c one recovers the equations 
(021 : as one should. Now, it is quite easy to present an 
example flow which the instability is triggered: take for 
instance a flow for which p is almost constant at infinity 
(almost means that it asymptotes to a constant value 
more rapidly than the sound velocity). Assume also that, 
near infinity, c = c\ + ^ . Then, we get that near infinity 
the effective potential behaves as 



2cic 2 k 2 



(52) 



For this to have a positive derivative, one requires C2 < 
(ci must be positive, as it is the asymptotic value of the 
sound velocity). We thus have one example of flow for 
which the instability is active. There are many others, of 
course, and there are also instances for which the system 
is stable. 



E. Absorption cross-sections 

The computation of absorption cross-sections may be 
handled analytically in the low frequency regime, to 
which I now turn. The computation of absorption cross- 
sections of different gravitational black holes has gained 
a special interest some years ago, since it was shown that 
string theory could reproduce these results (for some par- 
ticular geometries). I refer the reader to 18J for a intro- 
duction to the subject. Considering again our (2 + 1)- 
dimensional rotating acoustic black hole we shall now 
attempt at solving the wave equation in this geometry, 
in the limit of small u. The method we use here follows 
closely the work of Starobinsky and Churilov |5l| and 
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Unruh and others 66] . The wave equation reads |20| : 

2 



H 




v = f 

The computation will follow closely that in [63 . Chang- 
ing wavefunction to R = r~ 1 / 2 H we find that R satisfies 



Ad r (Ad r R)+ uj 2 r 2 - 2Bmu 



R = 0, 
(55) 

where A = r — 1/r. Let us solve this equation in the 
far-region, r >> 1. We then have, near infinity, 



Ad r (Ad r R) ~ r 2 d 2 r R + rd r R . 



(56) 



Assuming Blo << 1, the wave equation l|55|l in this region 
takes the form 

r 2 d 2 R + rd r R + (w V - m 2 ) R = . (57) 

Defining p = lot this takes the form 

p 2 d 2 p R + pd p R + (p 2 - m 2 ) R = , (58) 

which is a Bessel equation (see for example 0, ) , with 
the general solution 

R = aJ m (ojr) + pY m (ujr) . (59) 

Notice that m is an integer, and therefore J m and J_ m 
are not linearly independent. 

We will want later on to do a matching between the 
near region solution and the far-region solution, so let us 
investigate the near region behavior of this solution, or 
the limit tor — > 0. We find 63] 



R 



T[m + 1] 2 



7T Z 



(60) 



where ip is the digamma function. We now define the 
near-region as the range for which r — r + << — . Defining 



R = z a (l - zfT, 



(61) 
(62) 



then, in this region the solution representing ingoing 
waves at the horizon (which is the boundary condition 
one must impose) may be written as |H3] 



T = A\F[a, b, a + b - c + 1, 1 - z] 



where 



a = 1 + — — iw , 

m 

b = ivj , 

2 

c = 1 + m . 



(63) 

(64) 

(65) 
(66) 



Here w = , and F[., ., ., ] denotes the standard hy- 

pergeometric functions. Since c = m + 1 is a integer, 
one must be very careful in handling the hypergeometric 
function El El. 



When z — > 0. or r — > oo, we have 



R - Ar- m {i/j(a) + xjj(b) - ip(l + m) - + 



Ar m { '-""-'.TH , )■ (67) 
m\(a - m) m {b - m) m 



where A is some constant and {a) m stands for (a) m = 
a(a + \){a + 2)...{a + m— 1) , (a ) = 1. Matching the two 
solutions and l|H7jl we get 



(-l)"T[m] 2 (u;/2) 



-2m 



a _ 1 

n (a - m) m (b - m) m (ip(a) + ip(b) - ij)(m + 1) - tp(l)) 

(68) 

Now, (a — m)m — r[o— In] ' we ^ nus nave 



a _ 1 {-l) m T[m] 2 Y[a-m]T[b-m]{uj/2)- 2m 
/3 ~ ^r[a]r[6](^(a) + ^(6) - V(™ + 1) - ^(1)) 



(69) 



This is the final expression. Since a, /3 are related to 
the amplitude of in- and out-going waves at infinity, using 
(|69|l we can straightforwardly compute reflection coeffi- 
cients, absorption cross sections, etc. 



V. CONCLUSIONS 



Analogue black holes have proven to be a very valuable 
tool for the investigation of problems related to Hawking 
radiation. It is also possible that will yield valuable in- 
formation regarding classical phenomena involving black 
holes. We have shown here some aspects of classical phe- 
nomena involving acoustic black holes, that may prove 
useful for future experimental realization of these sys- 
tems. 
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